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ABSTRACT 

The present study addresses the subcritical 
aeroelastic response of swept wings, in various flight 
speed regimes, to arbitrary time-dependent external 
excitations. The methodology based on the concept of 
indicia! functions is carried out in time and frequency 
domains. As a result of this approach, the proper 
unsteady aerodynamic loads necessary to study the 
subcritical aeroelastic response of the open/closed loop 
aeroelastic systems, and of flutter instability, 
respectively are obtained. Validation of the aeroelastic 
model is provided, and applications to subcritical 
aeroelastic response to blast pressure signatures are 
illustrated. In this context, an original representation of 
the aeroelastic response in the phase - space is 
displayed, and pertinent conclusions on the implications 
of a number of selected parameters of the system are 
outlined. 

NOMENCLATURE 

a Dimensionless elastic axis position measured 

a 

from the midchord, positive aft 
c Chord length of wing, normal to the elastic 

axis, 2 b„ 

q Lift-curve slope 

h , h 0 Plunging displacement and its amplitude, 

respectively 
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/ , r a Mass moment of inertia per unit length of 

wing and the dimensionless radius of gyration, 
Jljmbl . respectively 

l Wing semi-span measured along the mid- 

chord line 

/ , m Dimensionless aerodynamic lift and 
moment, L a bJmU; and . 

respectively 

in , jj Wing mass per unit length and wing/air mass 
ratio, mjnpb] , respectively 
N Load Factor, 1 + h'/g 

p Peak reflected pressure in excess and its 

dimensionless value P m bJmUl , respectively 
r Shock pulse length factor 

s ,£ Laplace transform variable and operator, 
respectively 

c 7 Static unbalance about the elastic axis and its 

a ’Aa 

dimensionless counterpart, S y /mb * 
respectively 

/ , r Time variables and dimensionless time, 

U„t/b„ respectively 

1/ Freestream speed and its component normal to 

the elastic axis, respectively 
y Dimensionless free-stream speed ,U ft /b„CQ a 

Z Vertical displacement in z direction 

a,a 0 Twist angle about the pitch axis and its 

amplitude, respectively 
, £ a Structural damping ratio in plunging, 

c h /2mco h and in pitching, C(x /2I y a) a , 
respectively 

rj Dimensionless coordinate along the wing span, 

yfi 

A Swept angle (positive for swept back) 

£ Dimensionless plunge coordinate, wib n 
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X Dimensionless positive phase duration of the 

pulse, measured from the time of the arrival 
m Plunging-pitching frequency ratio ,(oJ(O a 

o,k Circular and reduced frequencies, cob ft !U n > 
respectively 

co h ,co a Uncoupled frequency in plunging, ficjm and 
pitching, ^jK a Jl y > respectively 

1. INTRODUCTION 

The modem post cold-war combat aircraft is likely to 
be exposed during its operational life to more severe., 
environmental conditions than in the past. This implies 
that, while being in flight, its structure should sustain 
pressure pulses due to fuel explosions, shock waves, 
sonic-boom, etc. It clearly results that the analysis of 
the aeroelastic response in the precritical range should 
be addressed in the various flight speed regimes i.e. 
from the incompressible to the hypersonic one. 

In this paper, using the time domain representation in 
conjunction with the concept of indicial function 1 “, the 
relevant unsteady aerodynamic loads necessary to 
approach the aeroelastic response of swept aircraft 
wings in various flight speed regimes, such as the 
incompressible, compressible subsonic, supersonic and 
hypersonic ones are obtained. Such a representation of 
unsteady aerodynamic loads is required towards, 
determination of the aeroelastic response to various 
pressure signatures 2 3 * 4 in the subcritical flight speed 
regime. Herein the idea of the modified strip theory 
initiated by Yates 5 has been further developed as to 
address the problem of the aeroelastic response of 
swept aircraft wings in various flight speed regimes. In 
this context, an original phase - space representation, 
provides full information about the behavior of the 
aeroelastic system. As a by-product of the response 
analysis, the conditions resulting in the occurrence of 
the flutter instability can be obtained. A validation of 
this procedure towards determination of the flutter 
instability boundary was done via the study of the 
aeroelastic eigenvalue problem, and an excellent 
agreement was reached. Herein the case of a swept 
wing (see Fig. 1), featuring the plunging and pitching 
degrees of freedom is considered. 

2. GENERAL THEORY 
2.1. Aerodynamic Loads 

An unified approach based upon the use of the 
aerodynamic indicial functions enabling one to 
determine the aerodynamic lift and moment in subsonic 
compressible and supersonic flight speed regimes is 
developed. The aerodynamic loads in the various flight 
speed regimes are obtained via the use of the 


appropriate indicial functions 67 and their application 
toward determination of the aeroelastic response and 
flutter instability boundary is addressed in the paper. 

2.1.1. Subsonic Compressible Flight Speed Regime 
A significant work toward developing appropriate 
analytical expressions for indicial functions in a 
subsonic compressible flow has been carried. 

In Refs. 3 and 4 it was shown that^ for subsonic 
compressible flight speed regime, the circulation 
around the airfoil is determined by a new set of indicial 
functions, different of the Wagner’s function. 




Fig. 1. a) Non-uniform swept wing; b) Airfoil 
section; c) 3-D view of a swept wing. 


In contrast to the incompressible case, the indicial 
functions in subsonic compressible flow are not 
analytic, except for limited instants of time. Following 
the formulation in Ref. 8, in our approach, new indicial 
functions, for the plunging and pitching degrees of 
freedom arc presented. Via the use of the four indicial 
functions defined for a two-dimensional lifting surface, 
the corresponding indicial lift and moment can be 
written as: 

L y ( t )= I C,„ pVlsWJbJtW' n a) 

M , (r) = C LaP VlSh' n (<p, M (r) + * (a + l)fc. (r)) , (lb) 

M r ) = P u - s A^ (rM(« + lk (r))’ (lc) 


b) 


I t* 
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M yil (T) = C La pU„S{2b)q 0 b 

(*«*,(*)+ i( fl + >X^c,(t)-0c« W)-7(« + 0 J ^„W).(ld) 
The expressions presented in Eqs. (1) have been 
adapted and extended as to approach the swept aircraft 
wing aeroelastic problems. Herein, the quantities in 
brackets identify the indicial functions for the 

compressible flight speed regimes; /i' is the plunging 
velocity and q {) is the indicial angular velocity; 

and {L q ,M yq ) correspond to the lift and 

moment due to plunging (in y direction) and pitching 
(around the elastic axis), respectively. The unsteady lift 
associated with the compressible flight speed regime 
can be described as: 



- 2 C La bpUi £ (r - r„ Kfo >/r 0 -(2) 

As concern the aerodynamic moment, M tl {z) this can 
be obtained from Eq. (2) by replacing and ^ by 
-2b<p cM and respectively. 

In order to have an unique formulation in both the 
incompressible and compressible flight speed regimes, 
the indicial functions appearing in the aerodynamic lift 
and moment, Eq. (2), have to be expressed in a form 
which is similar to that used in the incompressible flight 
speed range 8 , where the elastic axis is located at ah 
back from the mid-chord. To this end, their m odif i ed 
expressions have to be used. These are: 

0 < W-> 0 ,(r): (3a) 

(0 ^ <I>cm (0+ {(* + • V* ( r ) ; (3b) 

&,(O=> 0 «,( T )-{(« + l)M T ) : (3c) 

( r ) ^ (^ ) + 7 (« + ifytQ ( r )) ^3,J) 

-{(« + 

The validation of the new indicial function is provided 
in the frequency domain, where the unsteady 
aerodynamic coefficients arc determined, and the 
comparisons with those available in the classical 
literature 913 are provided in Fig. 2. The following 
representations for the unsteady aerodynamic lift is 
postulated: 

L a (k, t)= pUi^bi-L^ +olA (4a) 

b 


aerodynamic moment, A/„(r) this can be obtained from 
Eqs. (4) by replacing L t and //. by M, and A t > 
respectively. 
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where 

L } = iH x + // 4 , L 2 = iff 2 + 3 • (4b) 

Herein, the dimensionless unsteady aerodynamic 
coefficients //. have been introduced, and the reduced 

frequency k , has been included as to render the 
quantities in brackets dimensionless. As concern the 
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Fig. 2. Unsteady aerodynamic derivatives in: a-b) 
compressible subsonic and c-d) supersonic flight speed 
regimes. Comparisons with data from Ref. 9 through 
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Notice that there is a considerable influence of the 
Mach number on the indicial functions. In contra st to 
the behavior of the indicial function (Wagner’s 
function) in the incompressible flow speed that features 
at x - 0 an infinite value modeled as a Dirac pulse 5(r), 
in the compressible speed range a finite value for the„ 
indicial function at r = 0 is experienced. Moreover, in 
contrast to the compressible flight speed regime, where 
the aerodynamic lift and moment are obtainable., 
through four indicial aerodynamic functions, in the 
incompressible regime, a single function, namely the 
Wagner one is sufficient for determining the lift and 
moment. 

More recently, in Refs. 14 - 15, the concept of indicial 
function in subsonic compressible flow has been used, 
and an approximation and validation of indicial 
functions for any value of the Mach number in the 
compressible speed range was obtained. Keeping in 
mind the relationship between c(k) and for the 

incompressible flow-field (see Ref. 8), one can define, 
for the compressible flight speed regime, the analogous 
of Theodorsen’s function in terms of the corresponding 
indicial functions as: 

C c (k) = F c (k)+ iG c (k) = ik[ <j> c (r>'*Vr = ik& t (ik)\( 5a) 
C M (k)=M t (k)+iN c (k)=ik® M (; Ik ); (5b) 

c„ (*)■ K, (* ) + ‘G c , {k ) = ik<s> cq (i ik ); (5c) 

C rM „ (k)= M cq (k)+ iN" (k)= M<l (ik), (5d) 

where C c (k), C cM (k) are the compressible analogous 
of Theodorsen’s function for plunging for lift and 
moment and C c({ (k), C cMq {k ) are the compressible 

analogous of Theodorsen’s function for pitching for lift 
and moment, respectively. Also in this case, the four 
new complex functions, should collapse, in the 
incompressible flight speed regime, into a single one, 
namely the Theodorsen’s function. Paralleling the 
developments carried out for the incompressible flight 
speed regime, for the compressible one, the lift and 
moment per unit span for plunging (la/) and for 

pitching (l a/ ) about the leading edge, are expressed 
similarly as: 

L y (k)=C La pUlh'C c (k)-, (6a) 

M s (k)=2bC La pUlh'C M (*): C6b> 

L,,(0=2c la pt/>V l( ,(*); < 6c ) 

M „ (k ) = 2 C La pU„b 7 qC cUll (k). (6d) 

It clearly appears that, whereas from the response 
problem, considered in the time domain, the flutter 
boundary is obtainable in a numerical way, for the 
flutter analysis in a classical sense, the aerodynamic 
loads have to be converted from the Laplace space to 
the frequency domain, replacing s — > ik . 
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Fig. 3 Indicial lift and momentfunctions from plunging 
and pitching of a 2-DOF lifting surface. Comparisons 
with Refs. Xand 8. 
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displayed here. The same notations, as for the subsonic 
compressible have been used. The new supersonic 
indicial functions for selected Mach numbers are 
displayed and compared with those provided by Ref. 8 
in Fig. 3. Likewise Theodorsen’s function, also in this 
case, the newly "defined TndicTaT Tunctions can be 
expressed in terms of their real and imaginary parts. 
Some of these components are presented and validated 
TnTfgr47 In this context, the initial andTihaT values of 
the Wagner's function, expressed in the time domain, 
can be obtained by taking the limit of S $>(s) in the 
Laplace transformed space for r -> 0 and r — » °° , 
respectively. From the property of the initial value, in 
the incompressible flight speed regime, this implies: 
lim ^(r)= limsO(j) = <ft((T) = 0.5» 
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while, from the property of the final value: 

!im0(r)= lim$0(s)= 0(°°)= 1 * 

j ->0 

The first of the above expressions, namely the property 
of the initial value, is valid if <p{ r) does not contain an 
impulse $(j). In such a case the limit of the converted 

function in the Laplace transformed space does not 
exist. The second expression, occurring from the 
property of the final value, is applicable only in the case 
where 0(r) and 0'(r) are both Laplace transformable 

and if the function jO(s) has all its singularities in the 
left half of the s-plane. These properties are fulfilled by 
the Wagner’s function and by the new indicial function 
in the compressible flight speed regimes (see dotted 
fines In Figs. T.c and 4.d). 

2.1.3. Hypersonic Flight Speed Regime 
The supersonic flight speed regime is completed with 
that of the hypersonic flight speed regime, and in this 
context, piston theory aerodynamics is used (Fig. 5). 



analogous of Theodorsen’s function for the subsonic (a- 
b) and supersonic (c-d) flight speed regimes. 

2.1.2. Supersonic Flight Speed Regime 
The expressions of the indicial functions for two- 
dimensional lifting surfaces in supersonic flow are not 


Fig. 5 Aeroelastic response time-history in the 
hypersonic flight speed regime to blast loads, as 
presented in inset. Comparison of prediction of 
responses based on the present and piston theory 
aerodynamics (Ref. 16). 
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2.2. Sonic-Boom and Blast Pressure Signatures 

The blast and sonic-boom overpressure signature 
(referred to as the N-wave shock pulse) see Fig. 6 L can 
be represented mathematically in an unified way as 
(Fig. 6): 


rsv 


- EI^ + m^-S 1 ^--L t ,(y,t)=L l ,{y\t)’ <^ 9a ) 

dy 4 dr dr 

d 2 w T d 2 9 Ijr ^ ^ (%) 


67 If +Sy dt 1 


i y jpr-M a (y,t)=o- 




1--L 

t r J 






Herein, n{t) is the Heaviside step function; P m denotes 

the peak reflected pressure in excess to the ambient one 
(Ref. 3-4); t f) denotes the positive phase duration of the 
pulse measured from the time of impact of the 
structure; r denotes the shock pulse length factor. 
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For the cantilevered wing, the related boundary 

(8) conditions are: 

w ,(3T, f ) = M^ =0 (v,J =0’ and (10a) 

v=0 

3 2 w(y, /) _ 9 3 >v(y, t) __ dfl(y, /) _q, (10b) 

~W^~ 3? ” 3y 

Use of shapes functions, yields the aeroelastic 

governing equations in dimensionless form. In the right 
hand side member of these equations, L b (y,r) denotes 
externai 1 i me-dependent load act! ng on the rigid 
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Fig. 6 Blast and sonic-boom pressure signatures. 

For r = 1 the N-shaped pulse degenerates into a 
triangular pulse which corresponds to an explosive 
pulse, and for r = 2 a symmetric sonic-boom pulse is 
obtained. 

2.3 . Subcritical Aeroelastic Response of an Aircraft 
Wing Featuring Plunging-Pitching Coupled Motion 

Determination of subcritical aeroelastic response of 
swept lifting surfaces flying in a compressible flow 
field and exposed to time-dependent external pulses is 
useful in the design of wing structures and of the 
associated feedback control systems. The aeroelastic 
governing system of a swept metallic wing featuring 
plunging-twisting degrees of freedom and subjected to 
a blast pulse can be expressed as 8 : 


that can c orre s po n J ToTfFTast , sonic 
JfcHJom, gust or shock-wave pulse. Considering this load 
^tmiformly distributed in the cfiordwTse direction, ho 
moment contribution M b (y,r) is generated in Eq. (9). 
The aeroelastic governing equation can be converted in 

and solved for the 


the Laplace transformed space 
unknowns, | (=<£(£)) and «(=<£(«)); inverted back in 
time domain one obtain the plunging and pitching time- 
histories and the load factor time-history due to the 
sonic-boom pressure pulse. This inversion has been 
performed numerically via Mathematical routine 17 , and 
the results are depicted in the Figs. 7-10. 


3. RESULTS AND 1 


The graphs in Fig. 7-8 depict the dimensionless 


time-h istones for the plunging displacement % , the load 
factor N and the pitch 


a of a swept 
wing to a blast load in compressible subsonic and 
supersonic flow fields. In Fig. 7 the solid lines 
correspond to the coupled plunging-pitching motion of 
a 2 DOF model in the incompressible flight speed 
regime (M = 0) . 

The plunging-pitching coupling contributes to the 
reduction of the amplitude of the aeroelastic response. 
The load factor N has its maximum for t = 0, when the 
first impulse due the blast load occurs. 

In Fig. 9 a three-dimensional plot depicting the 
dimensionless plunging and pitching deflection time- 
histories £ of a swept aircraft wing (A = 15 e ) to blast 
pressure signatures vs. the normalized spanwise 
coordinate y and the dimensionless time r is supplied. 
The evolution of^The aeroelastic system can be 
graphically illustrated by examining its motion in the 
phase - space , rather than in the real space, and 
recognizing that die" trajectory depicted in this space 
represents the complete time history of the system. 
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Fig. 7 Aeroelastic response time-history of swept wing 
(a = 15°) in subsonic compressible flow field to blas t 
loads, as presented in inset 



Fig. 8 Aeroelastic response time-history of swept wing 
(a = 15°) in supersonic flow field to blast loads, as 
presented in inset. 
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Bunging Deflection 


V n =1.5 



Fig. 9 Three-dimensional plots depicting the 
dimensionless plunging and pitching a deflections 

time-histories of a swept aircraft wing (a = 15°) to blast 
pressure signature, vs. the normalized spanwise 
coordinate y and the dimensionless time r (M = 3). 

Corresponding to the condition resulting in the flutter 
instability (Fig. 10b and lib), that coincides with that 
obtained from the eigenvalue analysis, the trajectory of 
motion describes an orbit with constant amplitude 
{center). For V„ < Vf (Fig. 10a and I la) as time 
unfolds, a decay of the amplitude is experienced, which 
reflects the fact that in this case a subcritical response is 
involved ( stable focal point), while for V n > Vf (Fig. 
10c and lie) the response becomes unbounded, 
implying that the flutter instability has occurred 
{unstable focal point). 

To avoid its occurrence there are Iwo^ possibilities, 
namely, including a passive/active control methodology 
or acting on the sweep angle A. In the latter case the 
idea is to use the capabilities of the aircraft featurin g 
variable sweep angle (e.g. F-14 Tomcat) to reduce 
the severity of oscillations, and at the same time expand 
the flight envelope. 



Fig. 10 Phase-plane (§ ; §) portraits depicting the 
dimensionless plunging deflection time-history of a 
swept aircraft wing to a blast pulse, for selected values 
of the flight speed. 

The critical value of the flutter speed obtained from the 
eigenvalue analysis of the homogeneous system of 
equations, and of that extracted from the aeroelastic 
response coincide, being Vp = 2. On the basis of a 
n u mber of simulations, it can be co ncluded that the 
flutter predictions based on both methods are in an 
excellent agreement. 

The way to determine the flutter speed from response 
becomes clear fro m the phase plane portrait, in which, 
for a certain restricted range of the flight speed (in the 
vicinity of the flutter speed), a periodic response with 
constant amplitude is experienced. 
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4. CONCLUSIONS 




Fig. 11 Phase-space portraits (£;£ ; w) depicting the 
dimensionless plunging deflection time-history of a 
swept aircraft wing to blast pressure signature, for 
selected values of the flight speed. 


JThe versatility of the methodology presented here 
consisting of the approach of both the subcritical 
aeroelastic responses and of the flutter instability for 3- 
D swept aircraft wing in various flight speed regimes 
has been illustrated. The aeroelastic response has been 
represented in both the classical way, by displaying the 
time-histories of plunging - pitching motion and load 
factor, and in an original phase - space context, that 
provides full information about the behavior of the 
aeroelastic system. 

Applications assessing the versatility of this approach 
enabling one to treat both subcritical aeroelastic 
response a nd flutter instability were presented, and, as a 
continuation of this study, the treatment of actively 
controlled 3-D wing aeroelastic problems is 
contemplated to be achieved in the forthcoming 
developments. 
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